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Abstract 

The polynomial birth-death distribution (abbr. as PBD) onT = {0, 1, 2, ...} 
or X = {0, 1, 2, m} for some finite m introduced in Brown & Xia (2001) is 
the equilibrium distribution of the birth-death process with birth rates {oj} 
and death rates {/3i}, where > and [3i > are polynomial functions of 
i ^ I. The family includes Poisson, negative binomial, binomial and hyper- 
geometric distributions. In this paper, we give probabilistic proofs of various 
Stein's factors for the PBD approximation with ai = a and (3i = i + bi{i — 1) 
in terms of the Wasserstein distance. The paper complements the work of 
Brown & Xia (2001) and generalizes the work of Barbour & Xia (2006) where 
Poisson approximation (b = 0) in the Wasserstein distance is investigated. 
As an application, we establish an upper bound for the Wasserstein distance 
between the PBD and Poisson binomial distribution and show that the PBD 
approximation to the Poisson binomial distribution is much more precise than 
the approximation by the Poisson or shifted Poisson distributions. 
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1 Introduction and the main results. 



The cornerstone of the 'law of small numbers' is the Poisson limit theorem which 
says that the total number of independent (or weakly dependent) rare events follows 
approximately the Poisson distribution. Prohorov (1953), Hodges & Le Cam (1960) 
[also known as the Hodges-Le Cam theorem] and Chen (1975) can be regarded as 
the milestones for quantifying the Poisson limit theorem. However, Barbour and 
Hall (1984) proved that the accuracy of Poisson approximation in total variation is 
determined by the rarity of the events and its order does not improve when the sam- 
ple size increases. Inspired by the pioneering work in Poisson approximation, various 
attempts have been made to improve on the precision of approximation, leading to 
new forms of the 'laws of small numbers', e.g., binomial [Ehm (1991)], compound 
Poisson [Barbour, Chen & Loh (1992)], signed compound Poisson [Presman (1983), 
Kruopis (1986), Cekanavicius (1997)] and polynomial birth-death (abbr. as PBD) 
[Brown & Xia (2001)]. In particular, PBD approximation to the number of indepen- 
dent (or weakly dependent) rare events in total variation distance has been proved 
extremely accurate in Brown & Xia (2001), i.e. the errors of approximation decrease 
when the events become rarer and/or the sample size increases. On the other hand, 
all these approximations can also be considered in the context of the Wasserstein 
distance dw [see Barbour, Hoist & Janson (1992), p. 13]: for any two distributions 
V and Q on Z+ := {0,1,2,...}, 

' ' fdV- [ fdQ 



dw{V, Q) = sup 

where JF is the set of functions / : Z+ M, |/(x) — f{y)\ < |a; — ?/| for x,y E Z+. 
The metric dw characterizes the weak convergence and the convergence of the first 
absolute moments [see Shorack & Wellner (1986), pp. 64-65]. Poisson approximation 
with respect to the Wasserstein distance is well documented in Barbour, Hoist & 
Janson (1992), pp. 13-17 and Barbour & Xia (2006). In this note, we establish Stein's 
factors for PBD approximation in terms of the Wasserstein distance and the work 
generalizes the results in Barbour & Xia (2006). To demonstrate the significance 
of our results, we prove an upper bound for the Wasserstein distance between the 
PBD and Poisson binomial distribution. The bound is in the fashion of the Hodges- 
Le Cam theorem and it implies that the PBD approximation to the Poisson binomial 
distribution is much more precise than the approximation by the Poisson or shifted 
Poisson distributions. 

The PBD distribution introduced in Brown & Xia (2001) is the discrete distri- 
bution on X = Z+ or X = {0, 1, 2, m} for some finite m with probability function 




7rn = ¥^^\l + y^^^F^\ ,nGX, (1.1) 

where > and (3i>0 are polynomial functions of i G X. The framework for 
PBD approximation unifies Poisson approximation [Chen (1975)] adapted from 
Stein (1972) (in this case it is often appropriately called the Stein-Chen method), 
binomial approximation [Ehm (1991)], negative binomial approximation [Brown 
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and Philips (1999)] and some compound Poisson approximation [Barbour, Chen 
& Loh (1992)]. In fact, it is a routine exercise to check that a random variable X 
follows the distribution vr iff 

for every function g satisfying lE[/3x|5'(^)|] < oo. For any function / on T satisfying 
1/(0 < ^5 O'^s can recursively solve for the function gf such that 

Bgfit) := a,gfit + 1) - A^7/(0 = /(O " ^(/), (1-2) 

where 7r(/) = Xli^o /(O^*- Then for any random variable W, 

\EBgfiW)\ = \EfiW)-iTif)\ (1.3) 

gauges the difference between CW, i.e. the distribution of W, and vr in terms of 
the test function /. To estimate the Wasserstein distance between CW and vr, it 
is sufficient to calculate the supremum of the left hand side of (11.31) over the test 
functions / G JF. 

As in Brown & Xia (2001), we use PBD(a; 0,(3, 1) to stand for the distribution 
TT when ai = a and (3i = (3i + i{i — 1) with a, (3 > 0. Noting that the value of 
5f/(0) has no contribution to the equation (11.21) . we take 5'/(0) = 5^/(1) in this paper. 
In estimating the difference between the distribution of a random variable W and 
PBD(a; 0, (3, 1), one often needs to resort to the first and second order Palm distribu- 
tions of W [Kallenberg (1983), p. 103]. For the distribution of the sum of independent 
Bernoulli variables, it is shown in Brown & Xia (2001) that PBD(a; 0, /3, 1) offers 
very good quality of approximation. In fact, for random variables with well-behaved 
first and second order Palm distributions, PBD(a; 0, (3, 1) is a natural choice amongst 
many approximating distributions. 

Let II ■ II denote the supremum norm of the function ■ over its range. Depending 
on the choices of values of a and (3, we can often show that 

\^{ag{W + 1) - {f3W + W{W - l))g{W)]\ < e^gW + e^AgW + £2|| A^^^H 

for all functions g satisfying lE[iy^|5f(l^)|] < oo, where Ag{-) := g{- + 1) — g{-) and 
A'^g{-) := A(Ag{-)). It then follows from ([TS]) that 

dw{C{W),FBD{a;0,(3,l)) = sup |E/(1^) - PBD(«; 0, /3, 1)(/)| 

< eosup \\gf\\ +£:isup \\Agf\\ +£:2sup ||A^5(/||. 
/eJ^ /GJ^ /e.F 

Hence the major obstacle in applying Stein's method is on getting the right 
estimates of the solution g and its first two differences. The necessity of estimating 
the second difference of g here comes from the fact that the second order Palm 
distribution is used. We summarize the estimates in the following theorem, where 
?7i A r]2 := min{?7i, 772}. 
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Theorem 1.1 The solution gj to Stein's equation 

agf{i + 1) - [[3i + i{i - l)]gf{i) = f{i) - PBD(a; 0, /?, 1)(/) (1.4) 
with convention gf{0) = gf{l) satisfies 

II II ^V^^ 1 v/c. + (1-/3)74+ (l-/j)/2 ^^ 

k=l ^ a+2/3+2 

/GJ^ VP + l «(« + 2) + 2/3(a + /3 + 1)/ ' \ V a J 

2 13 
/g^" - a a + /5(l + /5 + a/2) - a ^ ^ 

Some care is needed when choosing the parameters in PBD distributions. For 
example, instead of using PBD(a; 0, /?, 1), one may prefer to use PBD(a; 0, 1, b) with 
ai = a and Pi = i + bi{i — 1) , i > 0, since, when 6 = 0, the distribution is re- 
duced to the Poisson distribution with mean a, denoted by Pn (a) . In other words, 
PBD(a; 0, 1, b) can be viewed as a generahzation of the Poisson distribution and its 
Stein equation becomes 

a~gf{z + + U{i - l)Ygf{i) = f{i) - PBD(a; 0, 1, b){f). (1.8) 

To relate to (II. 4p . one simply needs to define a = a/P, b = 1/P and cjf = Pgf. The 
following proposition is an immediate corollary of Theorem 11.11 

Corollary 1.2 The solution gf to Stein's equation U.8\) satisfies 

„ „ 1a \/ab+(b- 1)2/4+ (6- l)/2 , , 

sup||A^.|| < (^ + -, ^^^^T^ ,] /\(^ + -], {1.10) 

^ ~ \l + b a2 + 2a + 2 + 26(a + l); /\ Vv^ «/ 



2 b 

supllA^Ofll < - + - -. (1.11) 

/g3- ~ a (a + 1)6+1 + a/2 ^ ^ 

Remark 1.3 It is interesting to note that for Pn (a) approximation, i.e. 6 = 0, 
Barbour and Xia (2006) showed that the corresponding estimates of ^/ for Pn (a) 
are 

sup 11^/11 = 1, (1.12) 

n . „ 8 1.1437 

sup||A^;|| < lA— = < lA— ^, (1.13) 
/gj^ 3v2ea 

supllA^^^II < ^A-. (1.14) 
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It is straightforward to check that when 6 = 0, estimates (11 .Qp and (II. lip are reduced 
to (11.121) and (11.140 respectively, while (ll.lOp becomes f^^ff^, which is not the same 
as (11.131) . It is an interesting unsolved problem to estabhsh a bound of supjgjr || A^j|| 
for PBD(a; 0, 1, 6) which reduces to (I1.13P when 6 = 0. 

A prototypical example for applying the above estimates is to consider the Pois- 
son binomial distribution, i.e. the distribution of the sum W of independent Bernoulli 
random variables {Xi : 1 < i < n} with the distribution 

P(Xi = 1) = 1 - F{Xi = 0) =pi, l<i<n. 

Denote = J2^=iPiy = K/^i for r > 1 and A := Ai. Barbour and Xia (2006) 
proved that if A2 is an integer, then 

dw{C{W), Pn (A - A2) * 5aJ < 4(A - A2)"^A2= (1.15) 

where * denotes convolution and the point mass at A2. However, when A2 is not 
an integer, some minor adjustment is needed to make the argument work and the 
error bound will have to be increased by an amount reflecting the correction. The 
following theorem shows an impressive improvement on the quality of approximation 
by a PBD distribution since the bound in (I1.16P below is of order 63 + 6^ while the 
bound in (11.151) is of order 62 only. 

Theorem 1.4 With the setup in the preceding paragraph, set 

/5 = A^A^^ - 1 - 2A + 2A3A^\ a = (3X + X^ -X2. 
Provided A^Aj ^ — 1 — 2A > 0, we have 

dwiCW, PBD(«; 0, A 1)) < 3^3 + , , n ,n ■ (1-16) 

A — A2 — [i. + U2)U2 

Proof. The expansion of lEjBg{W) in Brown and Xia (2001), p. 1390 states that 

n 

i=l jj^i 
+ Y^P^PM+PJ)i^ -p^){l-pJ)n^'g{W'' + 1)], (1.17) 

where W' = W - Xi and W'^ = W - X^- Xj. Therefore, 

dw{C{W), PBD(a; 0, /5, 1)) = sup \Ef{W) - PBD(a; 0, (3, 1)(/) | 

= sup|E[%^(iy)]| 

< sup IIA^^jll Ip^P^i + ^p-Pj^(l -Pi-Pj) + ^PiPjiPi + Pj)^^- - Pi){l -Pj) 



i=l j^i j^i 



< (/3A3 + 2AA2)sup||A2(;^|| 

< 3{p\3 + 2AA2)/a, 



5 



where we applied (11.41) in the second equation, (I1.17P in the first inequahty and 
(11.71) in the last inequality. Since (3 > A^Aa ^ - 1 - 2A and A^ > A2, we have 
a= /5A + A2 - A2 > (A^Aa ^ - 1 - 2A)A + A^ - A2 = X^X2^ - A - A^ - A2 and a > f3X. 
Therefore, the above upper bound can be further estimated as follows: 

3/3A3 6AA2 ^ 3^ 6AA2 ^ 6g2A2 

a a - /3X + A^Aa'-A- A2-A2 ' A - A2 - (1 + ^2)^^2 ' ' 

Remark 1.5 If one uses (II. 6p and argues as in the proof of Theorem 3.1 of Brown 
and Xia (2001), the following estimate can be established: 

dwiCW, PBD(a; 0, P, 1)) < f ^ + f ^ + x O (^e,+VX9'A , 



where (Xfe = ^^YH=k+iPi^ Pi is the ith largest value of pi(l - pi), ^2(1 - P2), ■ ■ ■ , 
p„(l — Pn)- The bound here is not as good as the bound in (11.161) since the latter is 
of order 6*3 + 91 while a/A and x A/ \/X2 are usually large. 

Remark 1.6 Both PBD(a; 0, /9, 1) and shifted Poisson use two parameters, how- 
ever, the bound in (I1.16P is significantly better than (ll.lSp . In most applications, 
we can not expect A2 to be an integer which means some correction in shifted Pois- 
son approximation is unavoidable, while the PBD approximation does not have this 
inconvenience. 



2 The proofs. 

We follow the idea of Brown & Xia (2001) by setting gf{i) = hf{i) — hf{i — 1) for 
z > 1 so that 

Bgf{i) = ai[hf{i + 1) - hf{i)] + pi[hf{i - 1) - hf{{)] := Ahf{i), i > 0, (2.1) 

where A is the generator of birth-death Markov chains with birth rates {oi : i > 0} 
and death rates : i > 1}. In this setup. Stein's equation (II. 2p becomes 

and its solution can be exphcitly expressed as 

poo 

hfit) = - / E[/(Z,(t))-7r(/)]dt, (2.2) 

where {Zi{t) : t > 0} is the birth-death process with generator A and initial state 
i > 0. 

The advantage in estimating Stein's constants in Poisson approximation is that 
the couplings are relatively trivial and many computations are achievable [see Bar- 
bour and Xia (2006)]. For the PBD(a; 0,(3, 1) setting, we follow the main ideas in 
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Barbour and Xia (2006) but use immigration-death particle systems in {0, 1}^ to re- 
alize the couplings. Due to the generality of the model, our estimations are naturally 
more technical. 

Let Qk be the solution of (II. 2p for / = 1^, then Lemma 2.3 of Brown & Xia (2001) 
states that 

fi'fc(0 = , if fc < ^; gk{i) = n 5 It k>t, (2.3) 

aj_i7rj_i PiTTi 

where tt is defined in (11. ip . F{i) = Ylij<i'^j ^^"^ F{i) = Ylij>i'^j- 



2.1 The proof of dlS]). 

Noting that (7/(0) := (7/(1) and supjg^p \\gf\\ = supj>;^ supjgjr \gf{i)\, we will show in 
Lemma 12.11 that 

sup 1^/(01 = ^/i(2), 

where 

f^(k) = -k, ke Z+. 

We then show in Lemma [2.21 that (7/1 (z) is a decreasing function in i and reduce the 
estimate to 



sup 11^/11 =^/i(l) = - J^kTik, 



where the last equation comes from solving the Stein equation (ll.2p with i = 0. The 
estimate (/3 -|- ^+213+2 )'^ comes from (12. lip below. Now we concentrate on the other 
part of the upper bound in (ll.Sp . 

Noting that = a, Pk = Pk+k{k—l) and the balance equation a^Tik = Pk+i^^k+i 
holds for all k G , we have 

^kTTk < k'^TTk I = l^k{k - l)7rfc + ^ kiik 

k=l \k=l / \k=l k=l 

00 00 \ 1/2 / 00 00 \ 1/2 



,k=l k=l / \fe=l k=l 



^(/3fc - (3k)7rk + ^ /i;vrfc j = ( afc-iVTfc-i (1 - /5) ^ /ctt/, 

1 - /3) ^ fcvTfc 



\ 1/2 



A:=l 



hence 



fcvr, < (a + (1 - /3)V4)'/' + (1 - /5)/2. 

fc=i 



This implies 

/e^ « 
completing the proof of (II. 5p . 
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Lemma 2.1 For each i > 1, we have 

sup \gf{i)\ = gf^{i). 

Proof. Since the solution of (11. 2p does not change when / is shifted by a constant, 
it suffices to consider / G with f{i) = 0, that is, 

sup 15^/(01= sup 

From (12.31) . we know that gk{i) is positive for k < i and negative for k > i. Given 
that / G with f{i) = 0, it is easy to see that f{k) < \k — i\ and 

oo oo 

1^/(^)1 <^\i- k\gk{i) + - k\)gk{i), 

k:k<i k:k>i 

with the equahty attained at f{k) = i — k, or equivalently, at / = /i. ■ 
The following lemma summarizes the major properties of (7/^. 

Lemma 2.2 The functions gf^ and A^^fj^ are positive and decreasing in i >1, and 
Agf^ is increasing in i > 1. 

Proof. We will exploit a probabilistic interpretation oi gf-^ by constructing a coupling 
of {Zi{t) : t > 0} for all i G Z+ as follows [cf Barbour (1988) and Barbour & 
Brown (1992)]. Consider a particle system on the site space Z, the set of all integers. 
For each configuration ^ G {0, 1}^, set ^„ = 1 if position n is occupied by a particle, 
and if position n is vacant. Let {^(t) : t > 0} be a Markov process with state space 
{0, 1}^ and initial state ^n(O) = l{n>i}; which evolves as follows: 

1. particles immigrate to the system with rate a and a new immigrant takes the 
closest vacant site to the left of site 1; 

2. each existing particle suicides with rate /3; 

3. the particle at site i kills the particle at site j(> i) with rate 2; 

4. all above evolution rules apply independently. 

Fix i > 0. The evolution of J2n=-oo^n(t) does not depend on the particles at the 
sites to the right of i, i.e. {C,n{t) : t > 0,n > i}. It is straightforward to check that 
{Sn=-oo ^n{t) : t > 0} is a birth-death process with birth rate = a and death rate 
(3k = (3k + k{k — l), which means it has the same distribution as {Zi(t) : t > 0} := Zj. 
Take 

i 

m = ^«(*) (2-4) 



k=0 
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from now on. Let 

Ti := inf{t > : = 0} (2.5) 
be the first time that position i becomes vacant. 

Noting that gf{i) = hf{i) — hf{i — 1) with hf in (12. 2p . we have that, for any 
i > 1, 

POO 

Jo 

POO POO 

= E / ^i{t)dt = ETi = / P(Ti > t)dt, (2.6) 



which imphes that gfj^{i) > 0. Moreover, for any i >1, 

P(T, >t) = E[P(T, >t|Zi„i)] =Eexp|-^ (/3 + 2Z,_i(s)) rfsj , 

which is decreasing in i for any t since Zj_i(s) < Zi{s) for any i > 1 and s > 0. 
Hence gf^^i) decreases in i. 

Integration by parts ensures 

POO 

E(e-2(*^^»)|Z,_i) = / e-2(^^*)rf(-P(T, >s|Z,_i)) 

Jo 

= 1 - / 2e~^'F{Ti > s\Zi^i)ds. 
Jo 

Consequently, the increasing property of Agj^ follows from 

POO 

Agf,{i) = / (P(T,+i >t)- P(T, > t))dt 
Jo 

and 

P(Ti >t)- P(T,+i > t) = E (^e~ Io(P+^z,-i(-))ds (1 _ E(e-2(*^^')|Zi_i))) 

= E (^e-^'(^+2^'-(^»"'^^*2e-2T(T, > s|Z,_i)rfs^ 
= ^ ( Q- IoiP+^z,.^is))ds /"*2e-2^e-/o(/3+22i-iH)d«^5 



which is decreasing in i. Using the above equation again, we conclude that 
P(T, >t)- 2P(T,+i >t) + P(T,+2 > t) 

= E j\^-2s^- J^{f3+2Z,^^{u))du _ ^(^e~-2it^T,)~2{sAT,)^2,i^^)] ds^ . 

The decreasing property of A'^gf^ follows from the facts that 

POO 

A^gf,{i) = / (P(T,+2 >t)- 2P(T,+i >t) + P(T, > t))dt, 
Jo 
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is increasing in i and 

POO 

Jo 



-2{tAu)-2{sAu) 



) 



which is decreasing in i. Finally, A^gfj^ is positive because Agf/j is increasing. 



2.2 The proof of ([re]). 

Arguing as in subsection 12. 1^ we conclude that 



sup IIA^fjII = sup sup 



k=0 



sup sup 

i&+ /eJ^:/(i)=0 



fc=0 



Since Si = {Pi+i — Pi) — (aj+i — ai) = P + 2i is positive, Lemma 2.4 with condition 
(C4) in Brown & Xia (2001) states that Agk{i) is negative whenever k ^ i, which 
ensures that, for each j ^ T with f{i) = 0, 



\Agf{^)\<T.\'-''\^-^9k{^)) 



(2.7) 



k=0 

E (i - k)7lk 



F{i) F{i + l)\ >r^,, ., / F{i) 



a-Ki^i 



ani 



k:k<i ' k:k>i 

with the equality in (12. 7p reached when 

fik) = -\t-k\ = :Uk). 
However, since J2k£Z+ ^9kif) = 0, we have 

F{i) F{i + 1) 



TTfc 



k:k^i 



k:k<i 



ani_i 



aTTi 



E 

k:k>i 



F(.) 



F(2-r 



, F(i + 1) F(i-l) 1 1 

Q!^ Pj^ Oil Pi 



{21 



where the last inequality is due to Theorem 2.10 in Brown & Xia (2001). Now 
we compare the coefficients Ylk-k<ii'^ ~ ^)^« with Ylk-k<i''^k = F{i — 1), as well as 
Ylk-k>ii^ — i)7rk with F{i + 1). For the former, we use i — k < i. Apropos of the 
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latter, we apply the Caucliy-Bunyakovskii inequality to get 
^(A;-z)7rfe = F{i + 1)+ ^ {k - 1 - l)7rk 

k:k>i k:k>i+l 

1/2 / \ 1/2 



< 



\fc;A;>i+l / \k:k>i+l / 

/ \ 1/2 

\k:k>i+l / 



< F{t + l)+\ >^ P.n^j \lF{i + l) 

/ \ 1/2 

\k:k>i J 

= (1 + v^)F(z + 1). 
It now follows from (12.81) that 

F{i) F{i - 1) 



\^9nM < ^F{^-l)(^-^^i^)+{V^+l)F{^ + l) 



< max{z, + 1} {-Agk{i)) < max{z, + 1} f — A ) 



Hence, for i < y/a + 1, 



, , , , , a/o + 1 



and for i > ^/a + 1 , 



l^9/.>WI S ^ < < ^ < ^. 

which in turn ensure 

sup IIA^jll < + 

On the other hand, by (Q, O and (l23|) . 

A^7j,,(z) = E/ (|Z,+i(t)-z|-2|Z,(t)-z| + |Zi„i(t)-z|)dt 

^0 



/•oo 

E / i\Zi_i{t) - I + Ir.M + - 2|Z,_i(t) - % + 1t,m| + |^^-i(t) - «|)cit 

^0 

/•oo 

E / - z + 2| - 2|Z,_i(t) - z + 1| + - z|) lT„T.+i>trft 

+E / (|Z,„i(t) - 2 + 1|) (1t.m,t.+i<* - 

2 / P(Zi„i(t)=2-l,r„Ti+i>t)rft, 

r-oo 

+IE / (lz,_i(t)<i-l - lz,_i(t)>i) (lT,>t - lT,+i>t) dt. 

Jo 
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Now 

roc poo 

2 IP {Zi-i{t) =1-1, Ti, >t)dt<2 P {Ti, T,+i > t) dt 
Jo Jo 

1*00 i*oo 

Jo Jo 



00 



< 2 / e-2(/5+i)*dt 



+ 
and, by (12.61) and Lemma [2.2[ 

E / (lz,_i(t)<i-i - lz,_i(t)>i) (It.m - lT,+i>t) rft 
Jo 

^0 

< E / |lz._iW<i-i - lz,_i(t)>^| (P(T, > t|Zi_i) - P(T,+i > t|Zi„i))rft 

^0 

/"CX> 

< E / {T{Ti > t|Zi„i) - P(T,+i > t|Z,_i))rft 

Jo 

= Er,-Er,+i = -A(7^,(2)<-A(7^,(i). 

Solving the Stein equation (11.21) with i = and 1 gives 

/ N Ett , , 1 , , ^ , (a + /5)E7r — a 
gf,{l) = — , ^?^,(2) = - {Pi9f,{l) + Ett - 1) = ^ , (2.9) 

so — A(7jj(l) = ^ (1 — ^Ett). Collecting these estimates yields 

1 1 / /5E7r\ 1 a + 2 



(3 + 1 a\ a J- 13+1 a{a + 2) + 2/3{a + (3 + 1)' 
where the last inequality is from (12.111) below. ■ 

2.3 The proof of (1121). 

As in Brown & Xia (2001), we define r+ = inf{t : Zi{t) = i + 1} ioi i > and 
T~ = inf{t : Zi(t) = i — 1} for i > 1. For convenience, we set Tq" = 00, = Er^^ 
and e~ = Er~ for all i > 0. It has been proved in Brown & Xia (2001, p. 1378) that 

e+ = for z > 0, er = for z > 1 



and 



9,i^ = < '^^',1 (2-10) 



We begin the proof with two technical lemmas. 

Lemma 2.3 For i > 0, ~ '^'^t+i + ^ ^ ^.f^d for i > 1, e^_^2 ~ '^'^i+i + ^ 0- 
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Proof. For i > 0, 

F{t + 2) ^ Fjt + l) ^ Fjt) 

^ (A+2A+ii^(^ + 2) - 2A+i«F(^ + 1) + aTi^i)) 



^ / i+2 i+2 i+2 \ 

= — — A+2A+1 ^ TTfc - 2(3i+i ^ /^fcVTfc + ^ (3kf3k-iT^k I 

\ fc=0 A:=l A:=2 / 

= — ^ /3i+2/3i+ivro + (/3i+2/3i+i - 2/3i+i/3i)7ri + ^(/3i+2/3i+i - 2/3^+1/3^ + PkPk-i)'^k 

\ k=2 / 

Direct computation shows that /9i+2 > 2/9i and 

A+2A+1 ~ 2/5j+i/5fc + (3k/3k-i 
= A^i + (/^ + 2(^ + l))A+i - 2A+i/9fc + - + 2(fc - 1)) 
= (A+i - f3k? + + 2(^ + 1)) - (3k{(3 + 2{k - 1)) > 

for 2 < A; < i + 1, which imphes the claim for {ef : z > 0}. 

For 2 > 1, 

F(2 + 2) 2.^(^+1) , F(^) 



^ (A+iAF(z + 2) - 2AaF(z + 1) + ^^^(z)) 



^ 00 

= ^ V (A+lA - 2A/?fc + Pkh-l) TTk, 

* ^ A:=i+2 

SO it remains to show Pi+iPi — 2(3i(3k + PkPk~i > for A; > i + 2. Since 
I3k-i {k-l){(3 + k-2) f^_A(^ 1 



;3fc k{p + k-l) V V /3 + fc- 1, 

which is increasing in k, we can conclude that for k > i + 2, (3k-il3i+i > (3k(3i, which 
in turn gives 



A+iA - 2(3A + Mk^i > 2^ [3,+i(3i(3k(3k-i - 2(3i(3k > 0, 
completing the proof of the claim for {e~ : i > 1}. u 



Lemma 2.4 We have 

.2 



A%(1) < 



a/3 + 2{a + (3 + f3'^) ' 



Proof. Using (12.91) and the Stein equation (11.21) with « = 2, we have 
gjM = l(/?2^7^^(2) + E7r-2) 

= ^ ((2(/3 + + + a^) Ett - 2a{a + p + 1)) . 
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Using Lemma [2. 2 [ we have < (7/^(3) < gf-^(2), which ensures 

2a(« + /5 + l) <E.< ^(^ + 2/^ + 2) (2.11) 



2{p + l){a + p) + - - p{a + 2l3 + 2) + 2a 

Therefore, 

2 2{a + f3 + f3')En-2a{l + P) ^ 

^ fi'/iUJ - 5 S 



a3 - a(3 + 2{a + /3 + /3'^)' 

This completes the proof of the lemma. ■ 

The proof of ( HTTh . It follows from the definition of gf{0) = ^/(l), (EZD and 
Lemma 12.21 that 



|AV(0)| = |A^7;(1)| < A^;,,(l) = Agf.il) - 2A^o(l) < -2Ago{l) < ^ 



a 



where the last inequality is due to Theorem 2.10 of Brown & Xia (2001), so it 
remains to estimate A'^gf{i) for i > 1. In fact, we have from ( 12.10p that 

r -TTk (e+ 1 - 2e+ + e+i) , for A: > ^ + 2, 
^'^9k{i)=l TTi+i (6^+2 + 2e+ - e+ i) , foi k = i + 1, 
[ (er+2 - 2e,~+i + e,") , for A; < z - 1. 

Lemma 2.4 of Brown & Xia (2001) ensures that ef is increasing in i, so ej^_2 + — 
e^_-^ > 0, which, together with Lemma [2.31 implies 

A2^r,-^/-°' k < i-l or k = i + I, 

^9k{^)^ ifA:>z + 2. ^^'^^^ 

Replacing / by — / if necessary, it suffices to give an upper bound for A^gf{i). Again, 
we fix z > 0. Using / — f{i) instead of / if necessary, we may assume f{i) = 0, then 

= ^ f{k)A'gk{^) + E [f{k) - f{i + 1)]AV.(0 + /(^ + 1) E ^'9k{i 

k<i k>i+l k>i+l 

< - k)A^g,{{) + ^ [(z + 1) - fc]AV(0 + /(^ + 1) E ^'^'^(^)' 

fc<i k>i+l k>i+l 

with the equality reached when 

f{k) = {i- k)l{k<i} + [(z + 1) - A; + f{i + l)]l{fc>m} := fi^{k). 
If Efc>m^'^?^(0<0, then 

A^^^^jz) < - ^)A^^?fc(o + E [(^ + 1) - kW9k{^ - E ^'^^(^)' 

/c<i k>i+l k>i+l 
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with the equahty attained when fi3{k) = {i — A;)l{fc<j} + {i — k)l{^k>i+i} = i — k = 
fiik) +i, which is a shifted function of /i. On the other hand, if A'^g^.^i) > 0, 

then 

A'gjjt) < - k)A'g,{t) + ^ [(^ + 1) - k]A'g,{t) + ^'9k{i), 

k<i k>i+l k>i+l 

with the equahty achieved at fa{k) = {i — k)l{k<i} + {i + 2 — k)l{k>i+i} '■= faik). 
By Lemma [2.21 

AV.«<AVf.(l). 

On the other hand, since /j4 = 21{fc>j+i} + /i + i, we get from (12.121) and (12. 3p that 

oo 

A'gUi) = A'gf,{z) + 2 A'g,{^ 

k=i+l 

< A^gf,{l) + 2AViW = AVi(l) + 2(^?m(^ + 2) - 2g,+i{i + 1) + g,+,{ 
By Lemma [2. 3[ 

aiTi a-Ki-i ~ aTTi+i 

It follows that 

and the proof is completed by applying Lemma 12.41 ■ 
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